Combined holographic-mechanical optical tweezers: 
Construction, optimisation and calibration 

Richard D.L. Hanes, Matthew C. Jenkins, Stefan U. Egelhaaf 
Condensed Matter Physics Laboratory, 
Heinrich-Heine University, 40225 Diisseldorf, Germany 

July 21, 2009 



Abstract 

A spatial light modulator (SLM) and a pair of 
galvanometer-mounted mirrors (GMM) were com- 
bined into an optical tweezers set-up. This pro- 
vides great flexibility as the SLM creates an array 
of traps which can be moved smoothly and quickly 
with the GMM. To optimise performance, the ef- 
fect of the incidence angle on the SLM with respect 
to phase and intensity response was investigated. 
Although it is common to use the SLM at an in- 
cidence angle of 45 degrees, smaller angles give a 
full 27r phase shift and an output intensity which 
is less dependent on the magnitude of the phase 
shift. The traps were calibrated using an active 
oscillatory technique and a passive probability dis- 
tribution method. 



1 Introduction 

Since the invention of optical tweezers, they have 
had a major impact on the way nano-sized objects 
can be manipulated, especially in biology [2 El [3] 
and colloidal physics [H [S] . The first realisation was 
shown by Ashkin in 1970 where he demonstrated 
that spheres in a solvent of lower refractive index 
are subject to two forces. First, the gradient force 
draws particles into the centre of the laser beam 
due to the radial gradient in the beam intensity. 
Second, the scattering force accelerates the parti- 
cles in the direction of the beam due to radiation 
pressure. If the gradient force dominates the scat- 
tering force, which is achieved by tightly focusing 
the beam, a stable trap in three dimensions is cre- 
ated [7j , giving rise to optical tweezers [H [HI IB HH] ■ 



The trap can be moved within the sample by 
steering the laser beam using mirrors, which are 
typically mounted on galvanometers. Multiple- 
trap tweezers have been realised by splitting one 
beam into several traps in two or three dimen- 
sions using acousto-optic modulators [HI [12 [TH] 
or phase holograms [HI [El [El [HI [HI [IS] . Both 
of these techniques allow for motion of the traps, 
with their specific strengths and limitations [Hl20j. 
In addition to exerting a force on particles and 
thereby manipulating them, tweezers can also be 
used to measure forces in the pico- and femtonew- 
ton range [2T 1 [22 l [23] . 

This paper discusses a specific implementation 
of optical tweezers which combines a spatial light 
modulator (SLM) with galvanometer-mounted mir- 
rors (GMM). The SLM can generate complex ar- 
rays of traps in three dimensions, but the trap dy- 
namics are limited by the refresh rate of the de- 
vice. The GMM can translate the trap array in 
a fast and continuous way, increasing the flexibil- 
ity compared to an SLM alone. Furthermore, by 
combining both devices, multiple traps can be cal- 
ibrated simultaneously, allowing them to be used 
for force measurements. 



2 Optical tweezers set-up 

A schematic representation of the set-up is shown 
in fig [T] A diode-pumped solid-state laser with A = 
532 nm and a maximum continuous output power 
of 1.5 W (Laser Quantum Ventus 532-1500) is used. 
The polarisation angle is rotated by a half-wave 
plate (A/2), necessary for the SLM later in the opti- 
cal path. The beam is reflected from mirror Ml (all 
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Figure 1: Schematic diagram of the optical tweezers 
set-up combining a pair of galvanometer-mounted 
mirrors (GMM) and a spatial light modulator 
(SLM). (Abbreviations are defined in the main 
text.) 

mirrors from Thorlabs), then passes through a lOx 
beam expander BE (Newport NT55-578), which in- 
creases the beam diameter from 1.4 mm to 14 mm, 
the length of the longest edge of the SLM. It strikes 
the SLM (Holoeye LCR-2500) at an angle of inci- 
dence a = 22.5° and passes through the analyser 
(A) which, together with the half- wave plate (A/2), 
is necessary to optimise the device performance 
(section m. The light reflected by the SLM is then 
imaged onto the GMM by a telescope formed by 
lenses LI {fi—175 mm planoconvex, all lenses from 
Edmund Optics) and L2 (/2=100 mm planocon- 
vex) which also reduces the diameter from 14 mm 
to ^GMM — 8 mm. Between the two lenses the light 
is reflected vertically by mirror M2, to align the 
beam with the microscope epi-flourescence port. 
The GMM (Nuffield Technology Inc. Quantum 
Scan 30) consist of two galvanometer-mounted mir- 
rors with perpendicular axes of rotation for two- 
dimensional in-plane translation of the traps. From 
here, the beam traverses the second telescope (figH} 
formed by lenses L3 (/3=125 mm planoconvex) and 
L4 (/4=100 mm planoconvex) which reduces the 
beam to its final diameter of dobj = if 4/ fa) c^gmm 
= (6'GMM/6'obj) dcMM = 6.4 mm, to slightly over- 
fill the back aperture of the microscope objective. 
Combined imaging and trapping is facilitated by 
a dichroic mirror DM (transmission in the range 
525-540 nm is lower than 12 %, Chroma Technol- 




Figure 2: (Colour online) The rotation angle ^gmm 
of the galvanometer-mounted mirrors (GMM) re- 
sults in a tilt Oohj at the back aperture of the ob- 
jective. It is constrained by the geometry of the 
setup. Note that the GMM consists of two mirrors; 
only one is shown for clarity. 

ogy Corp. z532dcrb) which reflects the laser hght 
into a 1.4 NA lOOx Nikon objective. For imaging 
we use an inverted microscope (Nikon TE2000-U). 
The intense laser light is separated from the imag- 
ing light by a filter (notch filter, transmission in the 
range 525-538 nm is lower than lO"** %, Chroma 
Technology Corp. z532nf). The filter, the dichroic 
mirror DM and the lens L4 are housed within a 
home-built filter cube inside the microscope. The 
other optical components are mounted directly onto 
the optical table on which the microscope is also 
fixed. For capturing images a CMOS camera (Pix- 
elink PL-B742F) controlled by Lab View (National 
Instruments) is used. 

3 Galvanometers 

The galvanometer-mounted mirrors (GMM) are 
used to steer the beam into lens L3, part of the 
final telescope in the system (fig[2|). The telescope 
images the mirrors onto the back aperture of the 
objective, so that a rotation of the mirrors through 
an angle ^gmm corresponds to a rotation of the 
beam at the aperture through an angle ^obj, and 
an in-plane translation of the trap(s). Ideally both 
GMM would lie and rotate about a single point, 
but in practice this is not possible and they are 
thus shghtly separated (by 2.4 cm). This leads 
to a movement of the beam across their surfaces 
when they rotate and, as a result, the laser intensity 
through the objective varies. However, this varia- 
tion is less than 5% across the whole field of view, 
and all trap strength measurements were performed 
at the same position within the field of view. 
The diameter and focal length of lens L3 give rise 
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to a maximum useful mirror rotation angle ^q^m 
and hence a maximum translation distance within 
the sample. The diameter of lens L3 is 25.4 mm and 
the beam diameter at this point is doMM = 8 mm, 
giving a maximum displacement of the beam centre 
from the centre of lens L3 of (25.4 mm - 8 mm)/2 
= 8.7 mm. Furthermore, the distance between the 
central point between the two GMM and lens L3 is 
/3=125 mm resulting in a maximum rotation angle 
^'gmm = tan-i(8.7 mm/125 mm) = 4.0°. At the 
objective this corresponds to a maximum rotation 
^Sb]' = iU/h) ^GMM = 5.0° and an in-plane dis- 
placement of ±83 ^m. This displacement was de- 
termined experimentally and calibrated as follows: 
the filter was removed from the imaging path and 
the lowest available laser power used. This allowed 
us to observe the reflection of the beam from the 
interface between the cover slip and the solvent. 
The reflection was quantitatively tracked using a 
custom- written Lab View routine (section 15. ip . By 
moving a trap in the x and y directions the GMM 
voltage could be related to the lateral displacement 
in pixels (147 pixels/volt). No significant differ- 
ences were found for different positions or direc- 
tions within the field of view. The pixels were then 
related to absolute measures by comparison to the 
image of a Richardson test slide, which gave a value 
of 68.5 nm/pixel and hence a maximum displace- 
ment of ±83 fim. 

The GMM are controlled by a digital-to-analogue 
card DA (National Instruments PCI-6014) and 
home- written Lab View software. The output range 
from the DA card is -10 V to +10 V in steps 
of 5.8 mV, corresponding to a maximum rota- 
tion range of the mirrors of ±48° with incre- 
ment 0.03°. Since we are using a smaller range 
of angles |0gmm|<^gmm— ^.0°, we incorporated a 
lOx voltage divider between the computer and 
galvanometer control electronics (Quantum Drive 
3000) resulting in a maximum possible rotation of 
±4.8°>6lg^\i with increment 0.003°. Details of 
trap dynamics are discussed below f section 15. 3p . 

4 SLM characterisation 

To optimise the performance of the tweezers sys- 
tem the SLM was characterised in detail. Optimum 
performance entails a high relative intensity of the 
wanted first order projections with respect to the 



unwanted zero and higher order intensities, result- 
ing in the desired size and shape of each trap with- 
out the presence of ghost traps. An SLM working 
at optimum efficiency (ignoring the algorithm used 
to calculate the phase holograms [24]) will be able 
to modify the phase of the incoming light by 2tt, 
without significantly affecting its amplitude. Our 
SLM uses nematic liquid crystals to control the re- 
fractive index at each pixel, making it highly sen- 
sitive to the incoming polarisation orientation. It 
also causes a slight depolarisation which reduces 
the efficiency, and requires an analyser to recover a 
well defined polarisation. 

A Holoeye LCR-2500 SLM is used which gives 
a 2tt phase shift from 400 nm to 700 nm when il- 
luminated at a small angle of incidence (manufac- 
turer's specifications). For this device the effects 
of changing the incidence angle a on the output in- 
tensity and phase were investigated experimentally. 
The different a require interferometers with differ- 
ent geometries (fig The set-up used for zero 
angle incidence geometry represents a Michelson- 
Morley interferometer (fig [3^). The beam is 
first expanded by a beam expander BE and then 
split by a beam splitter cube BS (B.Halle TWK-20) 
between a mirror M and the SLM. The refiections 
are collected at the same BS, passed through an 
analyser and recorded on the camera. The incom- 
ing beam polarisation is rotated using a half wave 
plate (A/2). For incidence angles 3° < a < 22.5° 
we used another interferometer (figlSb). This set- 
up contains two beam splitter cubes, BSl to divide 
the incoming beam and BS2 to recombine them. 
A symmetric set-up with nearly equal path lengths 
is required due to the limited coherence length of 
our laser (about 1 cm). A similar set-up, but with 
rectangular geometry, a Mach-Zehnder interferom- 
eter '571 , is applied for an incidence angle a — 
45° (fig [St)* 

4.1 Output intensity 

Ideally, the SLM would introduce a phase shift 
without affecting the output intensity. Since this 
is in general not achieved for all conditions, the 
SLM was characterised to find its optimum config- 
uration. For different incidence angles a, we de- 
termined the output intensity / as a function of 
phase shift at the SLM for different polarisa- 
tions of the input ipi and output ipo beams. The 
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Figure 3: Interferometers used for determining the 
intensity and phase as a function of different inci- 
dence angles a. (a) a = 0°, (b) 3° < a < 22.5° and 
(c) a = 45°. For the intensity measurements, the 
beam directed towards the mirror M is covered and 
the camera used as a power meter. (Abbreviations 
are defined in the main text.) 



Figure 4: Output intensity / as a function of phase 
shift A'a0 introduced by the SLM for an incidence 
angle a — 22.5° and input and output polarisation 
ipi—lAQ° and (^o=0° respectively. 



mirror arm of the set-up was blocked using a beam 
dump, and the intensity / reflected from the SLM 
recorded with the camera. The intensity recorded 
across all pixels of the camera was averaged to im- 
prove statistics. It was ensured that when the laser 
was off the reading of the camera was zero and that 
saturation of the camera was avoided. The input, 
</Ji, and output, ipo, polarisations were selected us- 
ing the half- wave plate (A/2) and analyser. The po- 
larisations were defined relative to the laboratory 
vertical and increased clockwise in the direction of 
beam propagation. The phase shift A(/) introduced 
by the SLM was then increased in steps and the 
intensity / measured. An increment of 1 out of the 
256 levels that our SLM offers, i.e. < TVa,/, < 
256, was used. An example of the dependence of 
the intensity on the phase shift is shown in fig ID 

As a measure of the fiatness of the response, 
we use the ratio between minimum and maxi- 
mum intensity /min/^max, therefore values close 
to one indicate an output intensity / almost inde- 
pendent of the phase shift A0. (Note that some- 
times the output amplitude A is considered with 
/min//max = (^min/^max)^ EH]-) The depen- 
dence of /min / ^MAX on ipi and ipo is shown in fig [5] 
for different incidence angles a. We covered a range 
0° < ipi,o <180° in increments of 20° with the re- 
sults at (pi^o = 0° and (pi^o = 180° being identical 
within experimental uncertainties. 

The intensity maps (fig O show that there are 
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Figure 5: (Colour online) Ratio between minimum 
and maximum intensity /min/-^max (fig H]) as a 
function of input, Lp{, and output, ipo, polarisation 
for different incident angles a ~ 0°, 4°, 22.5° and 
45° (top to bottom). Plateau regions are observed 
at (v5i, tpo) ~ (160, 20) and (60, 100). 



clear plateaux with /min/-^max close to 1 for all 
values of a. However, a tweezer set-up with a — 
0° involves passing a BS twice which reduces the 
intensity transmitted to the objective to a quarter 
and a = 4° requires long beam paths increasing the 
sensitivity to vibrations. Both configurations were 
therefore rejected and only a —22.5° and 45° were 
investigated further. In addition to hardly affecting 
the output intensity (implying /min/-^max ~ 1), 
the SLM should also provide A0 > 2tt, with 
not significantly larger than 27r to maintain a large 
range of A^a^- This is discussed in the following 
section. 

4.2 Phase change 

We now relate A^a^ to the actual phase shift A4> 
(fig in]) and determine the maximum phase shift ac- 
cessible for the different parameters; ipi, (po and a. 
The same set-ups are used but with the beam dump 
removed from the mirror arm of the interferometer 
(fig [3]). The mirror M was aligned to produce ver- 
tical fringes with a fringe width of 64 pixels, which 
was shown to provide optimum conditions consider- 
ing processing errors and statistics due to the num- 
ber of detectable fringes on our camera (1280 x 1024 
pixels) [28j ■ The two halves (top and bottom) of the 
SLM are programmed separately, one being held at 
a constant N/^^=0, the other increased in steps of 
1. As the difference in A^a^ and hence the optical 
path length between the two halves of the SLM in- 
creases, the fringes detected on the two halves of 
the camera shift relative to each other. To recover 
the phase shift A0, we took cross-sections from the 
top and bottom of the image. By cross-correlating 
these two signals, and finding the first maximum, 
we recover the shift in pixels. This is converted to 
actual phase shift by multiplication by 27r and 
division by the period of the fringes in pixels, found 
from the Fourier transform of one of the signals. 

We investigated the maximum accessible range 
of the phase shift for incident angles a = 22.5° and 
a = 45° for the same range of input, (pi, and output, 
(fo, polarisations as above, i.e. 0° < pi^o < 180° 
in increments of 20° (fig [7]). Close to the opti- 
mum conditions, we used increments of 1° (data 
not shown). We found that, for a = 45°, there 
was no configuration which gave a full 27r phase 
shift with /min/-^max > 0.1. Therefore we use the 
a = 22.5° configuration which gives optimum re- 
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Figure 6: Actual phase shift A(/) as a function of 
N/^tf, for an incident angle a — 22.5° and input 
and output polarisation (pi=140° and iy9o=0° re- 
spectively. 

suits for ipi — 140° and (^o = 0°. These parameters 
result in a phase shift At/) = 2.07r and an intensity 
ratio Imin/Imax = 0.54. The actual phase shift 
A(f) as a function of N^^ for this configuration is 
shown in figure [S] The non-linearity is corrected 
for by the software. 

5 Trap calibration and use 

5.1 Active method: Oscillatory drag 
force 

We use an oscillatory calibration technique, based 
on the drag force method However, here, in- 
stead of the stage the SLM-generated trap is moved 
using the GMM, and a precise positional calibra- 
tion of the stage is thus not required. A triangu- 
lar wave signal with frequency = 1 Hz is sent 
to the galvanometers. The reflection of the laser 
from the cover glass-water interface was tracked 
using one of the pattern matching algorithms in 
Lab View. It convolves the image of the laser spot 
with a template image of the reflected laser, to 
determine potential locations of the laser beam. 
Subsequently, the locations are refined to sub-pixel 
accuracy. The size of the tracked object (in pix- 
els) determines the accuracy of the tracking as well 
as the maximum number of objects in the field 
of view. In our experiments the objects are typ- 
ically 60 x 60 pixels, resulting in an accuracy of 
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Figure 7: (Colour online.) Maximum range of the 
actual phase shift A(j) as a function of input, tpi, 
and output, <po, polarisation for different incident 
angles a = 22.5° and 45° (top, bottom). The cross- 
hatched areas represent configurations where the 
intensity of one of the fringe signals was too small 
to be measured, and therefore could not be de- 
termined. 
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15 nm, determined by tracking a reflected laser 
spot. With our camera (full frame 1280 x 1024 
pixels at 27 frames/s, or e.g. 100 x 100 pixels at 
500 frames/s) this results in a maximum of 360 ob- 
jects at a time resolution of 37 ms, or one object 
at a time resolution of 2 ms. The position of the 
laser x^it) and the voltage sent to the GMM V{t) 
were recorded and XL (t) oc sin~^(cos(27rj/(t-t-AtL))) 
and V{t) oc sin-i(cos(27rz^(i + Aty))) fitted to the 
data. The time delay Ato — A<l — Aiy provides 
the processing time of the hardware and software 
system; for a frame rate of 70 frames/s it is typi- 
cally 10 ms. We subsequently trap a particle with 
radius R and perform a similar procedure, but this 
time recording the particle position xp (t) and using 
xp{t) oc sin^^(cos(27ri/(i -I- Aip))) as a fit function 
to obtain the time delay Ati — Atp — Aty (fig [5]). 
The time lag At = Ati — Ato can be related to 
the trap stiffness k — QirriR/ At by balancing the 
drag force due to particle velocity u in a viscous 
medium (viscosity rf) F^i = dnriRv and the trap 
force -Ftrap — kAx — KvAt, assuming a harmonic 
potential. A typical example using polystyrene sul- 
fonate particles in water with R — 2.0 ± 0.1 /xm 
(Interfacial Dynamics Corporation), a laser power 
of 11.1 mW, at temperature T = 25± 1 °C (leading 
to 77 = (0.89 ± 0.02) X 10"^ Pa s) and a distance to 
the coverslip of about 15 /im is shown in figure [51 
The fit resulted in At = 20.8 ± 0.8 ms, and thus n 
= (1.6±0.1) X 10"'^ pN/nm and a maximum excur- 
sion of the particle from the trap center Axmax = 
134±5 nm. The assumption that the trap potential 
is harmonic up to Axmax was confirmed (section 
Oand fig©. 

The trap stiffness k — £,/At depends on the fric- 
tion coefficient 1^ = GnriR which above was cal- 
culated from ?7 and R. It can, however, also be 
determined independantly from the free diffusion 
coefficient of the particle in the bulk D — ksT/^ 
to give K = ksT/ DAt. D was determined from 
the mean squared displacement in two dimensions 
< >= ADt by recording the particle motion 
while diffusing. In order to track the particle accu- 
rately, it was kept close to the imaging plane and 
15 fim from the coverslip using the tweezers. The 
tweezers were programmed to trap and release the 
particle periodically with a frequency of 1 Hz which 
was realised by jumping the tweezers to a distant 
point using the GMM. To recover the free diffusion 
coefficient D, the data corresponding to the periods 
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Figure 8: (Colour online.) Voltage V sent to gal- 
vanometer {+) and particle position xp (o) together 
with fits as a function of time t. Polystyrene- 
sulfonate particles with a radius R ~ 2.0 /im in 
water at a distance to the coverslip of about 15 /im 
were subjected to an oscillatrory optical tweezer 
with a laser power of 11.1 mW at a temperature of 
T=25 ±1 °C. We determined a time lag between 
particle and trap position At = 20.8 ± 0.8 ms giv- 
ing a trap stiffness k = (1.6 ± 0.1) x 10~^ pN/nm. 



when the particle was influenced by the tweezers 
were discarded and the traces of the particle while 
it was freely diffusing shifted so that the data at 
the boundaries overlapped. To test this method we 
used the experimentally determined D to calculate 
the radius R = kBT/GirriD of the same particles 
as above (i? = 2.0 ± O.l/im manufacturers spec- 
iflcation's) making sure the selected particles are 
representative; very large or small particles were 
rejected based on their size in the recorded images. 
We obtain an average particle size of 2.04±0.06/im 
(which included the polydispersity of the sample) 
with an error bound of an individual particle ra- 
dius of about 2%. This compares very well with 
the manufacturer's specifications and indicates a 
successful determination of the friction coefficient 

5.2 Passive method: Probability dis- 
tribution 

A particle or several particles were held in sta- 
tionary traps and their positions determined us- 
ing video microscopy and particle tracking in two 
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dimensions. We record the position at 150 Hz 
for 300 s. This provides the probabihty distri- 
bution of the positions, p(r), whose bin size was 
chosen to coincide with the accuracy of the par- 
ticle tracking technique (fig [9|). A gaussian was 
fitted to p{r) and from the variance the trap 
stiffness k keT/cr^ determined [TT]. With the 
same particle and conditions as above we deter- 
mined 0-2 = (2.26 ±0.04) X 10"^^ m^ and thus k = 
(1.56±0.03) X 10"'^ pN/nm, which agrees within ex- 
perimental uncertainty with the value determined 
using the active method. 

5.3 Three dimensional particle ar- 
rays 

With our optimised SLM we can generate multiple 
traps (fjg fTO]) . The SLM allows us to trap particles 
in three dimensions. This represents an advantage 
over acousto-optic deflectors, with which it is pos- 
sible to trap particles in more than one plane [12] . 
but full control in three dimensions is not yet pos- 
sible. Our set-up furthermore provides improved 
flexibility in moving multiple particles due to the 
combination of the SLM with the GMM. This com- 
bination allows for fast and smooth movements of 
the trap array. The time taken for one galvanome- 
ter to move through 0.5° (equal to 10.5 fim in-plane 
translation) and settle to within 90% of its final po- 
sition is 0.65 ms (based on manufacturer's specifi- 
cations). With the SLM such a movement is not 
possible smoothly, due to its limited refresh rate of 
72 Hz (corresponding to 14 ms) for pre-calculated 
holograms, or, for real-time applications due to the 
limited processing speed of the host computer, a re- 
fresh rate of about 3 Hz (corresponding to 330 ms) . 

6 Conclusions 

We have constructed optical tweezers which com- 
bine galvanometer-mounted mirrors (GMM) and 
a spatial light modulator (SLM). We investigated 
the SLM with respect to its amplitude and phase 
shift and optimised the tweezers set-up accord- 
ingly. The optimum incidence angle for our SLM 
is considerably below the typically used 45°. With 
an optimised set-up we can create complex arrays 
of traps in three dimensions and move the array 
in a smooth and fast way using the GMM. The 
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Figure 9: (Colour online) Top: Two-dimensional 
probability distribution of particle positions p{r) 
with r — (x, y) within the trap. Contours of the 
data and fit at 5%, 20%, 40%, 60%, 80% and 99% 
of maximum (outside to centre) are shown. Bot- 
tom: Azimuthally averaged probability distribu- 
tion p{\r\) together with a gaussian fit with vari- 
ance (T^ = (2.26 ± 0.04) X 10"^^ m^, indicating 
a trap stiffness k = (1.56 ± 0.03) x 10"''' pN/nm. 
The inset shows the same data and fit in a log - 
square representation, and the vertical line shows 
that Axmax — 134 nm as used above is within the 
harmonic regime (see text for details). The same 
sample and experimental conditions were used as 
for fig [El 
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Figure 10: Image of six polystyrene sulfonate parti- 
cles with radius R = 2.0 fim, three trapped between 
the objective and the focal plane (e.g. top left) and 
three the other side of the focal plane (e.g. top 
right). 



traps can furthermore be used for force measure- 
ments. They were calibrated using two methods, a 
drag force technique and a probability distribution 
method, which yielded consistent trap strengths 
(within 3%). 
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